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Abstract 

We consider the possibility of observing continuous third-harmonic generation using a 
strongly driven, single-band one-dimensional metal. In the absence of scattering, the quan- 
tum efficiency of frequency tripling for such a system can be as high as 93%. Combining 
the Floquet quasi-energy spectrum with the Keldysh Green's function technique, we derive 
the semiclassical master equation for a one-dimensional band of strongly and rapidly driven 
electrons in the presence of weak scattering by phonons. The power absorbed from the driv- 
ing field is continuously dissipated by phonon modes, leading to a quasi-equilibrium in the 
electron distribution. We use the Kronig-Penney model with varying effective mass to estab- 
lish growth parameters of an InAs/InP nanowire near optimal for third harmonic generation at 
terahertz frequency range. 



Introduction 

When electrons in a crystal band are driven by an external time-independent electric field, they 
move periodically across the Brillouin zone, creating characteristic Bloch oscillations. 

[DEI Th e 

fre- 
quency of the oscillations, C0q = eEa/h, where a is the unit cell size, coincides with the energy 
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separation between neighboring states localized on a Wannier-Stark ladder.^ The effect has been 
observed for electrons/holes in semiconducting superlattices,^ for atoms trapped in a periodic 
optical potential, 9 and for light propagating in a periodic array of waveguides, with gradient of the 
temperature or of the refraction index working as an effective electric fieldJ^JS 

Combining the effects of a strong, time-periodic driving field, with the nonlinearity of the 
Bloch oscillations leads to higher harmonic generation of the driving frequency JEHU] This effect 
has recently been observed in bulk ZnO crystals strongly driven by a few-cycle pulsed infrared 
laser. 16 The application of the infrared field in short, 100-femtosecond, pulses was necessary to 
ensure that the absorbed energy could be transferred to the lattice and dissipated. 

In this work, we suggest that frequency multiplication due to periodically-driven Bloch oscil- 
lations could also be observed in a steady-state setting, e.g., a periodically modulated nanowire (or 
an array of such nanowires) continuously driven by high- amplitude terahertz radiation (see Fig. [I]). 
In the weak- scattering limit, the quantum efficiency of frequency tripling for such a system can be 
as high as 93%. 

For a nanowire in mechanical contact with an insulating, optically transparent substrate, a 
quasi-equilibrium electron distribution will be reached as the power absorbed from the driving 
field will be continuously dissipated into phonon modes. This distribution can be quite different 
from the initial, equilibrium Fermi distribution. In particular, at the driving field amplitude which 
is optimal for third harmonic generation, the distribution can be both broadened and inverted. 
The inversion of the distribution occurs once the driving field amplitude exceeds the dynamical 
localization threshold.* 17 * 18 ! 

In our analytic derivation, we combine the Floquet quasi-energy description with the Keldysh 
Green's function technique to obtain the semiclassical master equation for a one-dimensional band 
of strongly and rapidly driven electrons in the presence of weak scattering by phonons. We solve 
these equations numerically to find the electron distribution function for a cosine energy band 
at a given driving field frequency (fixed at £l/2n — 1 THz) and the field amplitude chosen to 
suppress the generation of the principal harmonic. This electron distribution is used as an input 
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(a) (c) 

Figure 1: (a) A nanowire made with alternating InAs/InP regions, (b), (c) Schematic of the third 
harmonic generation with a planar array of such nanowires. The driving field is ^-polarized so that 
the electric field E\ be parallel to the nanowires. The generated third harmonic will have the same 
polarization but propagate at a different angle. 

for calculating the time-dependent current and the intensity radiated at different harmonics of the 
driving field frequency. We use these results to find the optimal dimensions of a periodically 
modulated InAs/InP nanowire, which would yield the most efficient frequency tripling of 1 THz 
radiation. 



Theoretical approach 

We consider a single-band one-dimensional metallic wire driven by a harmonic electric field with 
the amplitude Eq and frequency £1, and coupled to substrate phonons, 

H = #0+#e-ph+#ph, (1) 

where the electron, electron-phonon, and phonon Hamiltonians are, respectively 

Ho = Y, e ( k+A (*)V k Ck, (2) 

k 

#e-ph = ^" V2 EM qjik cI + ^(fcq + fciq), (3) 
q,k 
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#ph = L^q^q. (4) 

q 

Here (<~[) is the annihilation (creation) operator for an electron with one-dimensional momentum 
fik and energy e{k). To apply our results to a periodically modulated nanowire, we assume a tight- 
binding model with the electronic spectrum, 

e(k) = -2Jcos(ka), (5) 

where / is the hopping matrix element and a is the period of the potential along the chain. The 
electric field is incorporated into the Hamiltonian through the vector potential A(t)=Ao sin £lt with 
Aq = eEo/h£l representing the vector potential of the driving field. Phonon annihilation (creation) 
operators b q and b\ are labeled with the three-dimensional wavevector q = (#|i,qj_) and co q is the 
phonon frequency (electron spin and phonon branch indices are suppressed). The factors M q ^ = 
a q fi(h/2co q ) 1 / 2 are the matrix elements for electron-phonon scattering. 

We ignore the effects of disorder or electron-electron interactions, and consider lattice phonons 
in thermal equilibrium at temperature ti/kBfi. We do not include directly the scattering by phonon 
modes of the nanowire, assuming that they are strongly hybridized with those of the substrate, with 
the corresponding effects incorporated in the matrix elements M q £. The electron-phonon coupling 
is considered to be weak, meaning that the phonon scattering time is long compared to the period 
T = 2n/£l of the driving field. 

Modified energy spectrum of the driven system 

The dynamics of the strongly-driven electrons with the Hamiltonian ([2]) is characterized by non- 
monotonous phases 

(fr(t) = J* dt'e(k+A{t'j). (6) 
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The phase accumulated over a period, %(t), can be expressed in terms of the average particle 
energy with the momentum hk, 

(e(k+A)) = t~ 1 dte(k+A(t)); (7) 

clearly, this energy can be also identified as the Floquet energy of a single-electron state. While 
Eq. ([7]) does not include the usual additive uncertainty m£l, this particular choice has the advantage 
that in the weak-field limit, Aq — » 0, (e(k+A)) recovers the zero-field spectrum e{k). 

The average energy ([7]) also coincides with that introduced in the theory of dynamical localiza- 
tion.'^^ Dynamical localization occurs when the effective band becomes flat, i.e., (e(k+A)) —> 0. 
The corresponding condition is most easily obtained in the special case of tight-binding model with 
the spectrum ([5]), 

(e(k+A)) = -2Jcos(ka), J = JJ (A a), (8) 

where Jo(z) is the zeroth order Bessel function. With the driving field amplitude increasing from 
zero the bandwidth is gradually reduced; it switches sign at the roots of the Bessel function, A§a — 
C$ n . The first time this happens corresponds to the electric field Eq = C$\MXjea, where £oi ~ 2.405. 

Frequency Multiplication with weak scattering 

We obtain the instantaneous current by averaging the canonical velocity operator dH/dA over the 
electron distribution function fa = (c^c^), 

i(t) = Cf{t) sinA(t)a + Sf(t) cos A(t)a, (9) 

where we assumed the tight-binding spectrum ([5]) and used the definitions 

C f (t) =^jj^ cos(fc) A, S f (t) =^J^ MQfk. (10) 
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In the limit of weak scattering, the distribution function fa is time- independent and always 
symmetric, = fa. Thus, Sf(t) = while Cf(t) = Cf is a time-independent pre-factor. The 
Fourier components of the current are obtained directly, 

i(t)=2Cf J m (Aoa)sm(m£lt), (11) 

ra=l,3,5,... 

where the summation is over the odd harmonics m. By choosing Aoa = £n w 3.8317, the first 
harmonic can be fully suppressed, which leaves the third harmonic dominant. The maximal value 
for the fraction of the energy emitted into the third harmonic (93.34%) is found in close vicinity of 
this amplitude, see Fig. [2} 



\Im? 
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Figure 2: Normalized magnitude squared of the Fourier harmonics of the instantaneous current, 
|/ m | 2 , for m = 1 (red dashed), m = 3 (black, solid), and m = 5 (blue, dotted) plotted as a function of 
the dimensionless amplitude of the vector potential of the driving field, see Eq. ([2]). The intensities 
\I m \ 2 correspond to the power emitted in the corresponding harmonics when multiple nanowires 
are used in a planar geometry, see Fig.[T](b),(c). 

Transition kinetics in a driven system 

We use the Keldysh non-equilibrium Green's function (GF) formalisirP2H22] along with a perturba- 
tion theory expansion with respect to the entire time-dependent electron Hamiltonian ([2]); the cor- 
responding evolution is solved exactly in terms of the phases ([6]). Previously, related approaches 
have been used, e.g., for describing ionization of atoms'^SEH anc [ the high-order harmonic genera- 
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tionElin t j ie f ultrashort laser pulses. Here, instead of solving the corresponding equations 
numerically, we take the limit of weak electron-phonon coupling and analytically derive the semi- 
classical master equation for electron distribution function averaged over the period of the driving 
field, see Eqs. ( fTT] ) and ( [T8] ). The same master equation can also be derived from the formalism 



by Konstantinov and Perel' 26 with the help of an appropriate resummation of the perturbation 
series.^ 

In the interaction representation with respect to the time-dependent Hamiltonian ([2]), the elec- 
tron operators acquire time-dependence e~ l(Pk ^Ck with quasiperiodic phases We separate these 
phases by defining the "lower-case" GFs 

8k(t2A) = e- i<p ^G k (t 2 A)e i ^\ (12) 



where the "upper-case" G^(^^i) is any of the conventional GFs introduced in the Keldysh formal- 
in QSE] j] iese ph ases introduce rapid oscillations in the self-energy, making the direct Wigner 
transformation difficult. We notice, however, that in the limit of weak electron-phonon coupling, 
the GFs ( [T2| ) are expected to change only weakly when both time arguments are incremented by 
the driving period T. This implies that in the following decomposition, 

gk(t2,h)=Y^gkAtjy~ imaT , (13) 

m 

t = t2 — t\is the "fast" time, while T = fa + h ) /2 is the "slow" time when it appears as an argument 
of thus defined Floquet components g^ m (t,T) of the GF. The Dyson equations for thus defined 
Keldysh gf m and retarded gf m GFs 20 have the form 

folb (14) 
8 mfi 8(t)+I? on , (15) 

where and are the collision integrals originating from the corresponding self-energy func- 



(id T +mQ)gl m (t,T) = 
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tions. The collision integrals being relatively small, both gf m and gf m are dominated by the m = 
components. 28 

To derive the semiclassical master equation, we write the equations for the m = components 
of the "lesser" g K and "greater" g > GFs,E21 perform the Wigner transformation replacing the fast 
time variable t by the frequency 0), and use a version of the Kadanoff-Baym approximation^ 

gf fl ((Q,T) = iA kfi ((D,T)f k (T), A kfi (co) « 5(©), (16) 

for the corresponding spectral function, A^co, 7") = ^[^o( C0 ' -0 ~~ §^ > o( a) ' ^)] = 23m < ?^ (o), T), 
where (T) is the non-equilibrium electron distribution function averaged over the period. This 
requires that the electron-phonon coupling be weak, and assumes that the electron spectrum renor- 
malization has been included in the Hamiltonian ([2]). 

The resulting master equation for weak electron-phonon interactions has the following standard 
f orm HH 

d „ , N f dk' 



dt 



fk(T) = J—{r k j[i-MT)]MT) 



-rvtfk>(T)[l-f k (T)]}, (17) 



where the transition rates are 



poo 

IV = 2^|S^(m)| 2 / dCQW^ico) 

x (n a + l)8(Aej$ - hot)) + n m 8(Aej$ + hco) 



(18) 



Here W] c ^((o) is the phonon spectral function (density of states weighted by the square of the 
coupling) for a given momentum q\\ = k' — k along the wire, see Eq. (|3|, n m = [exp(/3cu) — 1] _1 is 
the phonon distribution function, and the energy increment 



Aej^J = (e(k+A)) - (e(k'+A)) -mhQ, (19) 
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is the energy carried in or out by phonons, depending on its sign. Note that this energy includes 
m quanta of the driving field, emitted or absorbed, depending on the sign ofm = 0,±l,.... The 
matrix elements S^{m) are the Fourier expansion coefficients of the product of the two phase 
factors, 6?' 5 ^W-*'<%'M 9 where 8(p k (t) = <p k (t) —t(e(k+A)) is the periodic part of the phase. They 
satisfy the sum rule 

oo 

E \S k , k >(m)\ 2 = l. (20) 

m=—oo 

Clearly, the equilibrium Fermi distribution for fa is only obtained in the limit of small electric field 
amplitudes, such that Sk^{m) with m = gives the dominant contribution. 



Simulation results 

The results presented in this section have been obtained by numerically finding the stationary 
solution of the discretized version of the master equation ( fT7| ) with transition rates ( [T8] ). A simple 



model for the phonon spectral function, W^^co) = y 2 0(a) — s\k — k!\), was used, with the sound 
speed s — 5 x 10 3 m/s as appropriate for typical 3D acoustical phonons. Since we assume no other 
scattering mechanisms, the quasi-equilibrium distribution functions fa and other results do not 
depend on the magnitude of the electron-phonon coupling y 2 . 

We fix the phonon temperature at 4.2K, the lattice period a — 8.64 nm, the average electron fill- 
ing at 1/2 and choose the driving field frequency £lj2%= 10 12 Hz (energy h£l « 4. 14 meV). Also, 
the amplitude A$a — £n « 3.8317 is fixed, which corresponds to the point where the first harmonic 
generation is fully suppressed [see Fig.^. At this point the effective coupling is / = //o(Cn) ~ 
—0.4037, which creates an inverted and somewhat narrowed band. The effective bandwidth is 
smaller than h£l for J < 2.57 meV. 

In Fig. [3} we show the intensity I/3I 2 of the radiated third harmonic (in arbitrary units) as a 
function of the tight-binding hopping parameter J. The overall upward trend reflects the linear 
scaling of the current with J. The plot has a series of pronounced maxima and minima related to 
the structure of the distribution function fa, see Fig. [4j Indeed, at the first maximum of the radi- 
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ated intensity I/3I 2 , / = 2.7 meV, the distribution function has a well-defined minimum at k = 
and symmetric maxima dXk — ±71/ a [Fig. |4jb)]; notice the population inversion consistent with 
negative On the other hand, the distribution in Fig. |4^c) corresponding to the first minimum of 
radiated intensity, / = 4.5 meV, is much flatter. This flattening can be traced to a sharp increase 
of the transition rates connecting the regions of momentum space near k = and k — Tija. This is 
illustrated in Fig. [5} where transition rates between k — and k — %j a are shown. The correspond- 
ing phases 8(p K / a = — <5<po have only even harmonics m£l, m = 2, 4, . . ., and the threshold values of 
7 for different m correspond to sharp maxima of r 7r / a . 

2.5 
2.0 
1.5 
1.0 
0.5 

2 4 -r / An 6 8 

J (meV) 

Figure 3: Magnitude squared of the third harmonic of the instantaneous current (arbitrary units), 
plotted as a function of the tight-binding parameter 7, computed with Nk = 37,43,47,51 discrete 
momentum points as indicated in the caption. See text for other simulation parameters. The 
pronounced minima are caused by the flattening of the distribution function near thresholds of 
m-photon-assisted scattering between the vicinities of k = and k = ±71/ a, with m even. 

In Fig. [6j we show how the average power & radiated into the phonon modes scales with the 
tight-binding parameter J. While general dependence on J is monotonic, at J — 4.5 meV, where 
the third harmonic has a minimum, & changes slope. 

Proposed Experimental Design 

The simulation results in the previous section suggest that the optimal system for third harmonic 
generation would be a one-dimensional metallic conductor with an unrenormalized bandwidth 
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Figure 4: Solid lines: the stationary distribution functions obtained by solving discretized versions 
of Eqs. ( fT7| ), ( [T8] ) with Nk = 51 momentum points and the tight-binding parameters / as indicated. 
See text for other simulation parameters. Dashed lines: equilibrium Fermi distribution functions. 
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Figure 5: The transition rate (18) (arbitrary units) for scattering between the sites at k = — 71 /a and 
k = (black, solid) and the individual contributions from m-photon assisted processes as indicated. 
The vertical dashed line at / = 5. 1 meV indicates the threshold for the m = 2 transition, \4J\ = 2h£l; 
the peaks to the left and to the right of this point correspond to phonon emission and absorption, 
respectively. 



11 



4 1( ™ 6 8 

J (meV) 

Figure 6: Scaling of the average power (arbitrary units) dissipated into the phonon modes as a 
function of the tight-binding parameter The four curves on top of each other correspond to the 
same numbers of discrete momentum points Nk as in Fig. [3j 

close to 2.6 times the energy h£l of the driving field quanta (bandwidth of about llmeV for 
£1/2tc = 1 THz is needed), and a wide gap to reduce the absorption of the generated harmon- 
ics. One option to satisfy these requirements is to use modulated semiconductor nanowires. Here 
we estimate the growth parameters of an InAs/InP nanowire,^ which would have a near optimal 
band structure for generating the third harmonic of a 1 THz driving field. 

We calculate the band structure of the modulated nanowire modeling it as a stack of cylinders 
with isotropic (bulk) electron effective masses m^ nAs = 0.073m,, and mj nP = 0.021m e for the InAs 
and InP carriers respectively, as appropriate for the nanowire diameter we used. 31 We used the 
barrier height of Vb = 0.636 eV, found from the four-band model simulations, which is close to 
experimentally observed™! 0.6 eV. To ensure a relatively large gap, we chose the nanowire di- 
ameter d = 20 nm, and InAs well width w = 6.0 nm. Separating the radial and angular parts of the 
corresponding wave functions, we obtained a version of the Kronig-Penney model with effective 
mass modulation, and effective barrier dependent on the transverse momentum fiK n i. We plot the 
first few allowed energy bads as a function of InP barrier width b in Fig. [7} 

In particular, we conclude that an InAs/InP nanowire of diameter d — 20.0 nm, well width of 
w — 6.0 nm, and barrier width of b — 2.64 nm [Fig.[T](a)] would have the lowest band with a width 
of approximately 10.9 meV. The next band would be separated by a gap of 280 meV [Fig.|7J. These 
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Figure 7: Energies of the three lowest bands computed using the Kronig-Penney model with effec- 
tive mass modulation corresponding to an InAs/InP nanowire with diameter d = 20 nm, InAs well 
width w = 6.0 nm, plotted as a function of InP barrier width, b. The labels indicate the radial n and 
angular / quantum numbers of dimensional quantization. The dashed line at b = 2.64 nm gives a 
bandwidth 10.9 me V, or tight-binding parameter / = 2.7, corresponding to the first maximum of 
the third harmonic in Fig. [3j 

parameters are near optimal for third harmonic generation at £l/2n = 1 THz. 

One possible device design could involve depositing of a number of parallel modulated nanowires 
on a substrate, with an ^-polarized driving field incident on the surface at angle 6 so that the electric 
field of the wave be directed along the nanowires [Fig.[T](b),(c)]. Then both the reflected signal and 
the first harmonic are going to be propagating at the same reflection angle 0, while the propagation 
direction of the third harmonic can be found from the Snell's law, sin 6 = 3 sin0, which accounts 
for the wavelengths ratio. 

Discussion 

In this work we suggest a possibility that frequency multiplication due to periodically-driven 
Bloch oscillation may be possible in a quasistationary setting, with the help of a narrow-band one- 
dimensional conductor. A quasi-equilibrium electron distribution is possible because the energy 
absorbed from the driving field is continuously dissipated by the bulk phonons. 

For a periodically modulated InAs/InP nanowire with the period a = 8.64 nm, and the driv- 
ing field frequency £1/271 = 1 THz, the emission of the first harmonic is suppressed with the 
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dimensionless vector potential amplitude A^a pt 3.83, which gives the electric field amplitude 
Eq = hQ.Ao/e « 1.8x 10 6 V/m, corresponding to the energy flux of about 0.5 MWt/cm 2 . At this 
kind of power, many effects could lead to eventual run-away overheating of the system, e.g., direct 
absorption by the substrate, or even a relatively weak disorder scattering in the nanowire. We hope 
that a quasi-continuous operation would still be possible, with the driving field pulse duration of a 
few microseconds, as opposed to few picoseconds in the experiment. 16 
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